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Abstract 



We estimate the signal-to-noise ratio for two gravitational detectors inter- 
acting with a stochastic background of massive scalar waves. We find that 
the present experimental level of sensitivity could be already enough to detect 
a signal from a light but non-relativistic component of dark matter, even if 
the coupling is weak enough to exclude observable deviations from standard 
gravitational interactions, provided the mass is not too far from the sensitivity 
and overlapping band of the two detectors. 
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The sensitivity of present detectors to a stochastic background of rehc gravitational waves 
has been recently discussed in detail in many papers (see - [Q, for instance, and references 
therein). The sensitivity analysis has also been extended to include scalar waves ||^, and 
scalar stochastic backgrounds of massless (or massive, but light enough) scalar particles, 
interacting with gravitational strength with the detectors. At present, however, no analysis 
seems to be available on the possible response of the gravitational antennas to a scalar 
stochastic background of non-relativistic particles. 

The aim of this paper is to compute the signal-to- noise ratio (SNR) for a pair of gravi- 
tational antennas by taking into account the possible mass of the background particles, in 
order to discuss in some detail the possible effects of the non-relativistic branch of their 
spectrum. 

We shall consider a cosmic stochastic background of massive scalar waves, whose energy 
density is coupled to the total mass of the detector with gravitational strength (or weaker). 
We shall assume that the background is characterized by a spectral energy density fl{p) = 
d{p/ pc)/d In p, which we measure in units of critical density pc = SH^Mp/Sn, and which 
extends in momentum space from p = to p = pi (pi is a cut-off scale depending on 
the details of the production mechanism). As a function of the frequency f = E{p) = 
(m^ +p^)^/^, the spectrum fi(/), 

w).f|z^=(zys,(p) (1) 

thus extends over frequencies f > m, from / = m to / = /i = (m^ + pfY^'^ (note that 
we are using "unconventional" units in which h = 1, for a better comparison with the 
observable quantities used in the experimental analysis of gravitational antennas). We may 
thus distinguish three phenomenological possibilities. 

• m ^ /o, where /o is any frequency in the sensitivity band of the detector (tipically, 
if we are considering resonant masses and interferometers, /o ~ 10^ — 10'^ Hz). In 
this case we expect no signal, as the response to the background should be totally 
suppressed by the intrinsic noise of the detector. 

• m -C /q. In this case the detector, in its sensitivity band, responds to a relativistic 
frequency spectrum, and the SNR can be easily estimated by using the standard results. 
For a relativistic background of cosmological origin, however, the maximal amplitude 
allowed by nucleosynthesis is ~ 10~^, possibly suppressed by a factor -C 1 
(in the interaction with the antenna) to avoid scalar-induced, long-range violations of 
the equivalence principle (see for instance). We thus expect from such a scalar 
background a response not larger than from a background of relic gravitons, and then 
too weak for the sensitivity of present detectors. 
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• m ~ /q. In this case the mass is the frequency band of maximal sensitivity, and the 
detector can respond resonantly also to the non-relativistic part of the background 
(i.e. to the branch p < m of ^{p))- In the non-relativistic sector, on the other hand, 
the background amphtude is not constrained by the nucleosynthesis bound, because 
the non-relativistic energy density grows in time with respect to the relativistic one: 
it could be sub-dominant at the nucleosynthesis epoch, even if today has reached a 
near-to-critical amplitude ~ 1 (i.e., even if the massive background we are consid- 
ering represents today a significant fraction of the cosmological dark matter). In such 
case, it will be shown in this paper that the present sensitivity of the existing gravi- 
tational antennas could be enough to distinguish the physical signal from the intrinsic 
experimental noise. 

We will follow the standard approach (see [0], for instance) in which the outputs of two 
detectors, Sj(t), i = 1,2, are correlated over an integration time T, to define a signal: 

rT/2 

S= / dt dt'si{t)s2{t')Q{t-t'). (2) 

J-T/2 

Here Q{t) is a real "filter" function, determined so as to optimize the signal-to-noise ratio 
(SNR), defined by an ensemble average as: 

SNR = {S)/AS = {S) ((5^) - {SfY'^^ (3) 

The outputs Si{t) = hi{t) + ni{t) contain the physical strain induced by the cosmic back- 
ground, hi, and the intrinsic instrumental noise, nj. The two noises are supposed to be 
uncorrelated (i.e., statistically independent), {ni{t)n2(t')) = 0, and much larger in magni- 
tude than the physical strains hi. Also, the cosmic background is assumed to be isotropic, 
stationary and Gaussian, with (hi) = 0. It follows that: 

rT/2 

{S)= dt dt\hi{t)h2{t'))Q{t-t'). (4) 

J-T/2 

An explicit compuation of the strain, at this point, would require a specific model of 
the interaction between the scalar background and the detector. We will assume in this 
paper that the strain hi{t), like in the case of gravitational waves and Brans-Dicke scalars 
0], varies in time like the scalar fluctuation (j){xi,t) perturbing the detector (computed 
at the detector position x = Xi), and is proportional to the so-called "pattern function" 
Fi{n) = eab{n)D°/^ , where n is a unit vector specifying a direction on the two sphere. Cabin) 
is the polarization tensor of the scalar along n, and Df'' is the detector tensor, specifying the 
orientation of the arms of the i-th detector. 
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The field (f){x, t) may represent the scalar (i.e, zero helicity) component of the metric 
fluctuations generated by the scalar component of the background, as in or could even 
represent the background fleld itself, directly coupled to the detector through a "scalar 
charge" qi (for instance, a dilatonic charge), as discussed in To take into account this 
second possibility, we shall explicitly introduce the scalar charge in the strain, by setting 

hi{t) = qi(f){xi,t)eab{n)Df, (5) 

where = 1 for scalar metric fluctuations, and < 1 for long-range scalar flelds, phe- 
nomenologically constrained by the gravitational tests. The dimensionless parameter q^ 
represents the net scalar charge per unit of gravitational mass of the detector, and is in 
general composition-dependent 

To compute the average signal (^) we now expand the strain in momentum space, 

p=\p\, p/p = h, E{p) = f = {w? +p^f/'^, (6) 

{(Ph denotes the angular integral over the unit two-sphere), and we use the stochastic con- 
dition 

{ct)\p,h),ct){p',n')) = 5{p-p')5\n-h')^{p). (7) 

The isotropic function $(p) can be expressed in terms of the spectral energy density f2(p), 
defined by 

p = p,yt/lnpfi(p) = ^(|0p), (8) 
(Mp is the Planck mass) from which: 

By inserting the momentum expansion into eq. (^, and assuming, as usual, that the obser- 
vation time T is much larger than the typical time intervals t — t' for which Q 7^ 0, we finally 
obtain: 

{S) = q^q^-^ J J^^liP)QiP)m- (10) 

We have defined the overlap function 7(p) and the filter function Q{p), in momentum space, 
as follows: 
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IQn J 

Q(p) = [ dt'Qit - t')e2-^(p)(*-*'). (11) 



Note that the overlap function depends on the relative distance of the two gravitational an- 
tennas and on their particular geometric configuration. In the above equation, in particular, 
7(p) has been normahzed to the response of an interferometric detector to a scalar wave 

We need now to compute the variance AS"^ which, for uncorrelated noises, much larger 
than the physical strains, can be expressed as p[: 



/■r/2 

~ (5^) = / dtdt'drdr'{ni{t)n^{r)){n2it')n2ir'))Q{t - t')Q{T - t'). (12) 

J-TI2 



rr/2 

l-T/2 

It is convenient, in this context, to introduce the noise power spectrum in momentum space, 
Siip), defined by 



(n,(t)n,(r)) = 1 1 dp^,(p)e-2-^^(^)(*--). (13) 



Assuming, as before, that T is much larger than the typical correlation intervals t — t', t — t', 
and using eq. ( pTj) for Q{p), then yields 

T f dp 

P 

The optimal filtering is now determined by the choice (see 0| for details) 



AS' = ^ [ '^E{p)S,{p)S2{p)Q\p). (14) 
4 J V 



n( ^ \ 7(P)^(P) 

E^{p)Si{p)S2{p) 

where A is an arbitrary normalization constant. With such a choice we finally arrive, from 
eq. ([T0|) and (plSf), to the optimized signal-to- noise ratio: 



AS OTT^ 



pE^{p) 5*1 (p) 5*2 (p) 



1/2 

(16) 



It must be noted, at this point, that the functions Si{p) and ■y{p) appearing in the above 
equation are different, for a massive background, from the usual noise power spectrum Si{f), 
and overlap function 7(/), conventionally used in the experimental analysis of gravitational 
antennas. Indeed, S,j are defined as Fourier transforms of the frequency / = E{p), so that 
(see for instance eq. ([TB|)): 

J df 5,(/)e-2-/* = Jdp 5,(p)e-2-^(^')*, 

df 7(/)e-2"'^* = [dp 7(p)e-2"'^(f)*, (17) 
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from which 



Siip) = {df/dp)S,{f), 7(P) = {df/dpm). 



By introducing into eq. (|T6|) the known, experimentally meaningful variables Si, 7, and using 
/ = E{p) = [m^ we thus arrive at the final expression: 



SNR = qiq2 



57r2 



T 



(m2 +p2)5/2 s^(^^m^ + p^)S2{Vm^ + p') 



1/2 



(19) 



This equation represents the main result of this paper. For any given massive spectrum 
f2(p), and any pair of detectors with noise Si and overlap 7, the above equation determines 
the range of masses possibly compatible with a detectable signal {SNR > 1), as a function 
of their coupling to the detectors. 

For m = we have p = f, and we recover the standard relativistic result 0, modulo a 
different normalization of the overlap function. For m 7^ we shall assume, as discussed at 
the beginning of this paper, that the mass lies within the sensitivity and overlapping band of 
the two detectors, i.e. 7(771) 7^ 0, and Si{m) is near the experimental minimum. Also, let us 
assume that the non-relativistic branch of the spectrum, < p < m, is near to saturate the 
critical density bound < 1, and thus dominates the total energy density of the background 
(the contribution of the relativistic branch p > m, if present, is assumed to be negligible). 



To estimate the integral of eq. (|l^), in such case, we can thus integrate over the non- 
relativistic modes only. In that range, we will approximate Si and 7 with their constant 
values ai f = m. Assuming that the spectrum Q{p) avoids infrared divergences at p 
(like, for instance, a blue-tilted spectrum Q{p) ~ {p/piY, with 6 > 0), we define 

' dlnp n\p) = nl, (20) 







where fl^ < 1 is a constant, possibly not very far from unity, and we finally arrive at the 
estimate 



SNR ^ q,q,^^ 



m^Si{m)S2{jn) 



1/2 

(21) 



Following 0, the background can be detected, with a detection rate 7, and a false alarm 
rate a, if 



SNR > V2 (erfc"^2a - erfc"^27) . (22) 



For a first qualitative indication, let us consider the ideal case in which the two detectors 
are coincident and coaligned, i.e. ^ = 1, Si = S2 = S, qi = q2 = q, and the massive 
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stochastic background represents a dominant component of dark matter, i.e. Q^hioo ~ 1 
(where /iioo = -f^o/(100 km sec~^ Mpc"^) reflects the usual uncertainty in the present value 
of the Hubble parameter Hq). In such a case eq. (pi]), for an observation time T = 10^ sec, 
a detection rate 7 = 95%, a false alarm rate a = 10%, gives the condition: 

m^/^Sim) < ^10-3iHz3/2. (23) 
We will use here, for a particular explicit example, the analytical fit of the noise power 



spectrum of VIRGO, which in the range from 1 Hz to 10 kHz can be parametrized as [10 
S{f) = 10-^^sec 



3.46 X 10"^ ( — + 6.60 X 10"^ f f \ 



500 Hz/ V 500 Hz 

3.24 X 10^^ + 3.24 x 10 



/ f 



500 Hz 



(24) 



The intersection of this spectrum with the condition (^), in the plane |log S", logmj, is 
shown in Fig. 1 for three possible values of q^. The allowed mass window compatible with 
detection is strongly dependent on g^, and closes completely for < 10^'', at least at the 
level of the noise spectrum used for this example. We should then consider two possibilities. 

If the spectrum Vt{p) of eq. (|19|) refers to the spectrum of scalar metric fluctuations, 
induced on very small sub-horizon scales by an inhomogeneous, stochastic background of 
dark matter, then = 1 (since the detectors are geodesically coupled to metric fluctuations). 
In that case the detectable mass window extends over the full band from 1 Hz to 10 kHz, i.e 
from 10-15 to 10-" eV. 

If, on the contrary, scalar metric fluctuations are negligible on such small scales, and 
n(p) refers to the spectrum of the scalar background fleld itself, directly coupled to the 
detector through the scalar charge g, then this coupling is strongly suppressed in the mass 
range of Fig. 1, which corresponds to scalar interactions in the range of distance from 10^ 
to 10^" cm. Otherwise, such scalar fleld would induce long range corrections to the standard 
gravitational forces that would be detected in the precise tests of Newtonian gravity and of 



the equivalence principle (see |jTT| for a complete compilation of the bounds on the coupling, 
as a function of the range). 

Taking into account all possible bounds |jll|, it follows that, if the scalar coupling is uni- 
versal (i.e. the induced scalar force is composition-independent), then the maximal allowed 
charge q^ is around 10~^ from 1 to 10 Hz, and this upper bound grows proportionally to 
the mass (on a logarithmic scale) from 10 to 10^ Hz. Composition-dependent couplings are 
instead more strongly constrained by Eotvos-like experiments: the maximal allowed value of 
q^ scales like in the previous case, approximately, but the bounds are one order of magnitude 
stronger. 
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FIG. 1. Tie boid curve corresponds to the noise power spectrum of VIRGO given in eq. (p^j. 
The thin, dashed Unes represent the minimal sensitivity required for the detection of the background 



at the 90% confidence level (i.e., SNR ~ 2.5), according to eq. (23). The mass window compatible 
with detection corresponds to the range of frequency for which S is below a given dashed line. 



By inserting into the condition (^) the gravitational bounds on we are led to the 
situation illustrated in Fig. 2. A scalar background of nearly critical density, non- universally 
coupled to macroscopic matter, turns out to be only marginally compatible with detection 
(at least, in the example illustrated in this paper), since the line of maximal is just on 
the wedge of the noise spectrum (0). If the coupling is instead universal (for instance, like 
in the dilaton model discussed in ^j), but the scalar is not exactly massless, then there is a 
mass window open to detection, from 10~^^ to 10~^^ eV. 

It seems appropriate to recall, at this point, that it is not impossible to produce a cosmic 
background of light, non-relativisic particles that saturates today the critical energy bound. 



as shown by explicit examples of spectra obtained in a string cosmology context [0. Such 
particles, typical of string cosmology, are in general very weakly coupled to the total mass 
of the detector (like the dilatons, if they are long range, and the charge of the antenna 
is composition-dependent), or even completely decoupled (like the axions, since the total 
axionic charge is zero for a macroscopic, unpolarized antenna). Nevertheless, it is important 
to stress that they could generate a spectrum of scalar metric fluctuations, gravitationally 
coupled to the detector, which follows the same non-relativistic behaviour of the original 
spectrum. We know, for instance, that in cosmological models based on the low-energy 
string effective action, the variable representing the dilaton fluctuations exactly coincides 
with the scalar part of the metric fluctuations (at least in an appropriate gauge [0), and 
that the associated spectra also coincide. 

In view of the above discussion, the results illustrated in Fig. 1 and Fig. 2 suggest a new 
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log^^Q (m/Hz) 

FIG. 2. The noise spectrum of Fig. 1 is compared with the maximal values of (as a func- 
tion of mass) allowed hy gravitational tests, in two cases: composition-dependent and composi- 
tion-independent scalar interactions. The thin dashed lines corresponds, from left to right, to 

= 10~^, 10^'', 10^^, 10^^, 10~^. The region compatible with a detectable signal is above the 
noise spectrum and below the bounds given by the gravitational experiments. 

possible application of gravitational antennas, which seems to be interesting. Already at the 
present level of sensitivity, the gravitational detectors could be able to explore the possible 
presence of a light, massive component of dark matter, in a mass range that corresponds to 
their sensitivity band, in spite of the fact that such a massive background could be directly 
coupled to the total mass of the detector with a charge much weaker than gravitational, or 
only indirectly coupled, through the induced spectrum of scalar metric fluctuations. 
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